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Abstract 

It is shown that viscoelastic wave dispersion and attenuation in a viscoelas- 
tic medium with a completely monotonic relaxation modulus is completely 
characterized by the phase speed and the dispersion-attenuation spectral 
measure. The dispersion and attenuation functions are expressed in terms 
of a single dispersion-attenuation spectral measure. An alternative expres- 
sion of the mutual dependence of the dispersion and attenuation functions, 
known as the Kramers- Kronig dispersion relation, is also derived from the 
theory. The minimum phase aspect of the filters involved in the Green's 
function is another consequence of the theory. Explicit integral expressions 
for the attenuation and dispersion functions are obtained for a few analytical 
relaxation models. 
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Notation 



z complex conjugate of z 

]a, 6] {x I a < X < 6} 
{a} one point set 

= P for t > and otherwise 

characteristic function of the segment [a, &] 

/(a;) ~a 5(2;) asymptotic equivalence lim^j^a [/(a;)/.9(a;)] = 1 

/(x) = Oa[5(a;)] < lim^^„ [/(x)/g(x)] < 00 

/(x) = Oa[g(x)] lim^^a [/(x)/g(a;)] = 

D" /(x) derivative D" /(x) = d" /(x)/dx" 

/'(a;) fix) - D/(x) 
u J := du/dt 

M, K+, C real, positive real, complex numbers R+ := {x G M | x > 0} 

C+ right half complex plane C+ = {z G C | | argz| < 7t/2} 

upper/lower half complex plane := {z G C | ±3z > 0} 

Lf — f Laplace transform of / f{y) — fi^) exp(— yx) dx 

Laplace transform of a Radon measure {Lfi){x) := ^, e~^^ /i(dt/). 



1. Introduction. 

The assumption that the relaxation modulus of a viscoelastic medium 
is a completely monotonia function has many implications for dispersion, 
attenuation and the Kramers-Kronig (K-K) dispersion relations. These as- 
pects are very important for acoustics, in particular for the ultrasound ap- 
plications. Ultrasound attenuation and mechanical tests are used as com- 
plementary methods of investigating molecular relaxation in polymers and 



soft matter [33| . Interconversion of ultrasonic data and mechanical test data 
is important for improving the reliability of the results. In materials with 
dielectric properties molecular relaxation can be investigated by studying 
dielectric loss. Dielectric loss in polar dielectrics exhibits similar features 
(alpha and beta peaks) to the mechanical loss modulus while the dielectric 
permittivity, representing gradual polarization of the material in an electric 
field, has all the properties of mechanical creep 

In seismological applications the minimum-phase aspect of viscoelastic 
wave propagation is relevant for deconvolution. 

The assumption that the relaxation modulus is a completely monotonia 
function has its roots in the interpolation of experimental data in terms of 
Prony sums with positive coefficients [l],[2|- By Bernstein's theorem complete 
monotonicity of the relaxation modulus is equivalent to the statement that 
the relaxation spectral measure is non-negative. A Prony sum with positive 
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coefficients is a completely monotonic function with a spectral measure sup- 
ported by a finite set of points. Several attempts have been made to justify 
complete monotonicity of the relaxation modulus by an assumption about 
dependence of stress on loops in the past history jsf, by an assumed relation 
between time dependence of strain and stress [3] or by other arguments jsj. 
The theory of viscoelasticity based on the assumption of complete mono- 
tonicity of the relaxation modulus is very attractive because it leads to a 
fairly complete characterization of the creep compliance, complex modulus, 
attenuation and dispersion as well as the possible anisotropic properties of 
the medium [6|. A weak point of the theory is an instability of the prop- 
erty of complete monotonicity: in an arbitrary neighborhood of a completely 
monotonic function in the space of bounded continuous functions there are 
functions which are not completely monotonic. The assumption of complete 
monotonicity has therefore to be considered as an a priori restriction on the 
data and on the interpolating functions. Parameter estimation based on ex- 
perimental data has to be performed in the class of completely monotonic 
functions. This is however implicit in routine modeling of experimental data 
for stress relaxation in terms of Prony sums for general viscoelastic media ^ 



or in terms of Cole-Cole 25, 26 , Havriliak-Negami 31-331 and Kohlrausch- 
Williams- Watts |35| relaxation moduli for more specific classes of viscoelastic 
materials. The same relaxation functions are used in modeling dielectric re- 
laxation in the same materials if they exhibit dielectric properties 0, 0] . 

Completely monotonic relaxation moduli and the associated complex 
moduli are determined by a relaxation spectral measure, which represents 
the weight of participating Debye relaxation mechanisms (Sec. [3]). It is shown 
below that the complex wave number function, the attenuation function and 
the dispersion function of a medium with a completely monotonic relaxation 
modulus are parameterized by the dispersion-attenuation spectral measure 
(Sec. ED. 

The dispersion function and the attenuation function satisfy a dispersion 
relation in parametric form because both functions are expressed in terms of 
a single spectral measure. This dispersion relation implies the K-K disper- 
sion relations with two subtractions (Sec. [TU]) . The K-K dispersion relations 
indicate that the complex wave number function is the Laplace transform 
of a causal distribution. This distribution turns out to be a second-order 
distributional derivative of a causal function defined in terms of the Laplace 
transform of the dispersion-attenuation spectral measure. 

The K-K dispersion relations are closely related to the minimum phase 
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property of viscoelastic Green's functions (Sec. [T^ . This property is relevant 
for deconvolution of seismic signals. 

The dispersion-attenuation spectral measure can often be explicitly cal- 
culated for a given analytic complex modulus by analytic continuation to the 
entire complex plane cut along the negative real semi-axis and calculating 
the jump of the analytic continuation at the branch cut. This procedure is 
demonstrated for the Cole-Cole, Havriliak-Negami and Cole-Davidson relax- 
ation models (Sec. [T2|) . The dispersion-attenuation spectral measure provides 
an efficient tool for numerical determination of the dispersion and attenuation 
functions. 

2. Mathematical preliminaries. 

Definition 2.1. A function f defined on the open positive real semi-axis M_|_ 
is said to be completely monotonic (CM) if it has derivatives of arbitrary 
high order and 



A CM function can have a singularity at 0. A bounded CM function / 
has a limit at and its domain of definition can therefore be extended by 
continuity to the closed positive real semi-axis. A CM function is locally 
integrable if and only if it is integrable over [0, 1]. 

Theorem 2.2. (Bernstein's Theorem) Every CM function f has an integral 
representation 



where ^ is a positive Radon measure (that is, a locally finite measure), such 
that 



for some e > 0. 

The Radon measure fi is uniquely defined by f . 

Definition 2.3. Wl is the set of positive Radon measures ^ satisfying the 
inequality 



(-l)"D"/(a;) > 



for a; > and n = 0, 1, 2, . . . 




t > 



(1) 




[0,oo[ 




(2) 
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Theorem 2.4. A real-valued CM function is locally integrable if and only if 
the Radon measure fi in eq. ([T]) belongs to Wl. 

The proof is given in Appendix I D[ Locally integrable CM functions will 
be denoted by the abbreviation LICM and the set of LICM functions will be 
denoted by £. 

Remark 1. Inequality is satisfied if and only Jjg /i(dr) < oo and 
Jj-^ ^jyu(dr)/r < oo hold simultaneously. This statement follows immediately 
from the inequalities 1/2 < 1/(1 + r) < 1 and l/(2r) < 1/(1 + r) < 1/r 
holding on [0,1] and on [l,oo[, respectively. By a similar argument is 
equivalent to the inequality 



for an arbitrary number a > 0. 

Definition 2.5. A function f defined on the closed positive real semi-axis is 
said to be a Bernstein function if f > and D/ is completely monotonic. 

Since a Bernstein function / is non-decreasing and non-negative, the limit 
limj._j.o+ f{x) always exists and is finite. 

A superposition of Bernstein functions with positive coefficients is a Bern- 
stein function. The set 53 of Bernstein functions is closed under pointwise 
convergence [I^ . Consequently an integral of a family of Bernstein functions 
with a non-negative weight is a Bernstein function. 

A function g is a. Bernstein function if and only if it has the form 



where a > and /i is a LICM function. Indeed, g' = h is a. CM function and 
g > 0, hence g is a Bernstein function. On the other hand, if is a Bernstein 
function then its derivative g' is a CM and it is locally integrable. Hence 
eq. ([3]) holds with h := g', a = g{0). In particular the primitive function 
1 — e~^ of the CM function e~^ is a Bernstein function. 

Definition 2.6. A function f : IR+ M. is a complete Bernstein function 

(CBF) if there is a Bernstein function g such that f{x) = g{x). 




]0,oo[ 




(3) 
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The set of complete Bernstein functions will be denoted by the symbol 

n. 

Let g be the Bernstein function ([3]). 
The Laplace transform of g has the form g{x) = a/x + h{x)/x, where a > 
and h is a LICM function. By Bernstein's theorem h{y) = J^^^^e'^^ fi{dr), 
where fi G VJl. Hence 

J[0,oo[ + ' 

If / is a CBF and /(x) = x^^(x), then 

f{x) = ax + x[ f^i^''") ^ x>0 
J[o,oo[ x + r 

or, equivalently, 

f(x) = ax + b + x[ t^^^ a; > (4) 

where fi G 971, a > and b := /i({0}) > 0. Any function with an integral 
representation of the form (jlj) with a,b > and ^ & Wl is obviously a CBF. 
In particular, the function x/{x + a) is a CBF if a > 0. 

We now note that x/ (x+a) is a Bernstein function because it is an integral 
of Bernstein functions x — )■ 1 — e"''^, r > 0, with a positive weight e~"^: 



X 

a 



x + a 



/ e-"'' (1 - e"™) dr 
Jo 



Eq. (ll]) now implies that every CBF function is an integral of Bernstein 
functions with a positive weight, hence it is a Bernstein function. 

The Bernstein function f{x) := 1 — e~^ is however not a CBF. Indeed, 
suppose the contrary. Then f{x) has the integral representation (jlj) with b = 
/(O) = and a = lim^^oo = 0. Hence f{x)/x = J^^^^^^^x{dr) / {x + r). 
But f{x)/x = e~^y xio,i]{y) therefore the characteristic function 

X[o,i] (x) = J]Q^[ e"^^ /u(dy) is a smooth CM function. This conclusion is false, 
hence 1 — is not a CBF. Two examples of complete Bernstein functions 
relevant for attenuation and dispersion are and (1 + x)" — 1, < a < 1. 

Many other examples of CBFs can be found in [10]. 
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Let g{z) denote the analytic continuation of g{x) to the complex plane 
cut along the negative real axis. Eq. (jl]) implies that '^g{z) > in C'^. By 
the Pick-Nevanlinna theorem (jlo[ Theorem 6.7) every non-negative contin- 
uous function g{x) on R+ which has an analytic continuation with the above 
property is a CBF. This criterion allows identifying some complex analytic 
functions as analytic continuations of CBF, in particular with < a < 1 
and ln(l + z). 



We shall also need the following theorems [10|. They follow easily from 
the Pick-Nevanlinna theorem. The first one is an immediate consequence of 
the Pick-Nevanlinna theorem. 

Theorem 2.7. If f is a CBF and < a < 1 , then /(■)" is a CBF. 

Theorem 2.8. A function f ^ is a CBF if and only if the function x/ f{x) 
IS a CBF. 

Proof. Let g{x) := x/f{x). 

If / is a CBF then it has an analytic continuation f{z) to C"*". The 
analytic continuation of / has an integral representation of the form (j4]). It 
follows that f{z)/z swaps C"*" and C~. Consequently its inverse z/ f{z) maps 

and into itself and is non- negative on R_|_. Hence g{x) = is a 

CBF. 

Conversely, /(x) = x/g(x). Consequently, if g{x) is a CBF, then /(x) is 
a CBF. □ 

For simplicity we shall henceforth apply the term CBF to the analytic 
continuation of a CBF as well. 



3. Viscoelastic media with non-negative relaxation spectrum. 

We shall henceforth assume that the relaxation modulus G is LICM joj. 
Since the function G is non-increasing and non- negative, it has a limit Goo '■= 
limt_!.oo G{t) > 0. The limit at can be infinite even for physically realistic 
models such as the Rouse theory of dilute polymer solutions [18|]. 

According to eq. 

G{t) = I e-^V(dr), t > 

J[0,oo[ 
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where /i G DJl. The relaxation modulus is thus a linear superposition of De- 
bye relaxation functions with non-negative weights. The relaxation spectral 
measure /i is positive. 

On account of eq. (jl]) the function 

Q{p):=pG{p)=p [ ^ (5) 

i[o,oo[ P + r 

is a CBF. It follows from the general theory of the Laplace transform (ill, [l^ 
that (5(0) = Goo > and, additionally, \imp^ooQ{p) = Gq := limt^Q^ G(t) if 
the relaxation modulus is bounded. 



4. One- and three-dimensional viscoelastic Green's functions. 

Viscoelastic Green's functions are solutions of the problem 

pu,u = G{t)*V^Ut, u(0,x) = 0, Mt(0,x)=5(x) (6) 

The viscoelastic Green's function in a three-dimensional space can be 
expressed in terms of the one- dimensional Green's function: 



«(3)(t,x) = --^f««(t,r) 
znr or 



(7) 



r=|x| 



Indeed, we have 

1 1 rico+£ poo 1 

^' > 2n2mJ_,^^, Q{p) + K{pf 

hence 

u^'\t,x) = — / e^' . f . . e-(^)Ndj9 (8) 
where the complex wave number function k{p) is defined by the equation 

<p)=p'^'p/QipY^' (9) 

The square root is defined in such a way that ^n{—iu) > for real u. 
In the three-dimensional case 

1 1 /-ioo+e /• ik-x 



M^^^(t,x) = — — — / eP*-f-dp / dgA; 



(271)3 27ti7_io,+, g(p) -y^3 p + ^(p)2 

1 1 /-ioo+e /-oo /-TT pifelxl cos(i9) 

/ eP*-^dp / A;2dA; / sin(?9) dt9- — 
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Hence, by closing the contour over k in the upper half of the complex fc-plane 



_-| rico+e poo iA:|x| 

^' ' (271)3 |x| g(p) ^7_^p + ^(p)2 



100+£ 



87t2i|x| Q{p) 2nrdr 



r=|x| 



5. The complex wave number function. 



Assume that under constant strain the stress does not relax to 0: Goo : = 
hm^^oo G{t) > 0. 

By Theorems 12 . 71 and 12 . 81 the complex wave number function k is a CBF. 
Since Q{0) = Goo > 0, k(0) = 0. Hence k has an integral representation (jl]) 
without the constant term: 

k{p) = Bp + (3{p) (10) 

where B > and 

W:=p/ ^ (11) 

J]o,oo[ P + r 

and p e DJl. The function /3 will be called the dispersion- attenuation function. 
The dispersion- attenuation spectral measure v determines the attenuation, 
dispersion and wavefront behavior of viscoelastic wave motion. 

We have just shown that the complex wave number function of a vis- 
coelastic medium with a LICM relaxation modulus is a CBF. The converse 
is not true. Indeed, Q{j)) = plp/i^ip)]"^ is a square of the CBF function 
p^^"^ p/ k{p). a square of a CBF need not however be a CBF, for example p"^^^ 
is a CBF but its square is not even a Bernstein function. 

Theorem 5.1. If^p>0 then ^f3{p) > 0. 

Proof. If > then 

m?) = [ t'lf Kdr) > (12) 

i]o,oo[ \p + r\^ 

□ 

Theorem 5.2. 3?/3(p) is a non- decreasing function of \p\ in the closed right 
half plane C-j.. 
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Proof. The integrand of eq. fll2p has the form f{x) := (x^ + ax)/(x^ + ax + c), 
where x = \p\, a := rcos(arg(p)) > 0, c = > 0. The derivative of / is 
non-negative. □ 

The derivative 



J]o,oo[ (p + ry 

exists for real p because the integrand is bounded by z/(dr)/(j9 + r). We are 
however interested in growth of ^/3{p) as p tends to infinity in the closed 
right complex half-plane. 

Theorem 5.3. //z/([0,oo[) = oo then ^(3{p) —)■ oo as |p| oo in the right 
half of the complex p-plane. 

Proof. The integrand of (|T2l) tends to 1 as \p\ — )■ oo. Setting x = l/\p\, 
a := cos(arg(p)) the integrand is transformed into the function (1 + ax) /(I + 
2ax + x'^), which has a non-positive derivative. Hence the integrand of f lT2|) is 
monotonically increasing to 1. If z/([0, oo[) = oo, then by the Lebesgue-Fatou 
lemma [l^ 5^/3 (p) increases to infinity as |p| — t- oo, p G C+. □ 

Theorem 5.4. f3{p)/p tends to as \p\ oo in the right half-plane —ti/2 < 
argp < 7z/2 uniformly with respect to argp. 

Proof For -n/2 < (f < n/2 

inf |i?e^'^ + rP = inf + + 2Rr cos(ip)] = R"^ + 

(pe[-7r/2,7r/2] (pe[-7t/2,7t/2] 



hence 



sup 

<fie[-n/2,n/2] 



13 {Re'^) 



< 



z/(dr) 



[o,oo[ vR^ 



r 



2 



But for R > l/^/3 and r > 1 the inequality 1/VR^ + r^ < V^/VTTS^ < 
1/(1 -|- r) holds in view of the inequality (1 + r)^ < 1 + 3r^ for r > 1, 
whi le for R > ^3 and r < 1, 1/VR'^ + < l/VS + r^ < 1/(1 + r). Thus 
l/VR"^ + < 1/(1 + r) for > 73 and r > 0. In view of eq. (E]) 



z/(dr) 



lim 



by the Lebesgue Dominated Convergence Theorem. Hence \(3{p)/p\ tends to 
in C+ uniformly with respect to argp G [—n/2, n/2]. □ 
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Corollary 5.5. If Gq < oo then B = (p/Gg)^^^ eg. (fTUi) . otherwise 5 = 0. 
Proo/. 

Recall that limp_5.oo = Go < oo. 
Theorem 15.41 implies that 



□ 



The phase function in the Green's function ([8]) has the form —iut — 
ia;)r = — iu (t— Br)— (3 {—iu) r with (3 {—iu) = Ooo[i^] and | exp(— r)| < 
1. It will be shown that Green's function vanishes for t<BriiB>0 and 
thus Coo := can be identified as the wavefront speed. The parameter B, 
B > 0, will henceforth be replaced by Coo varying in the range < Coo < oo. 

Eq. (ITOl) implies that the function k can be analytically continued to the 
complex plane cut along the negative real axis, i.e. to the principal Riemann 
sheet p G C, —n < argp < n. The jump of the complex analytic function 
on the branch cut ] — oo, 0] can be easily calculated. Note that the function 
X{p) := K,{p)/p assumes the following boundary values on both sides of the 
branch cut: 

^ ^ Jo r-i?±iO+' 

and 

{z ± iO+)"^ = vp^"^ =F 'm6{z) 
where vp z^^ denotes the principal value of z~^ and 6 denotes the Dirac 



delta [ij]. Assume for a while that z/(dr) = h{r) dr, where the density h is a 



smooth function. From the above identities follows the equation 

f{R) := [A (i?e'") - A (i?e-'")] /(2i7r) = h{R) (13) 
The boundary values of analytic functions are in general distributions 



15l |. Hence the jump function / is a distribution. Since / is a non- negative 
distribution, it is a measure. More generally the measure z/ of a segment 
]u,w] is given by the formula 

1 /""' 

iy(\u,w]) = - Q\ (i?e'") dR (14) 

J u 

(Theorem I G.3p . 
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6. A necessary and sufficient condition for complete monotonicity 
of the relaxation modulus. 

Theorem 6.1. The function K,{p) of a viscoelastic material satisfies the con- 
dition i^{pY /p E O. if the relaxation modulus is completely monotonic. 

Proof. Eq. ^ implies that Q{p) = pp^ / i^{pY. In a viscoelastic medium with 
a completely monotonic relaxation modulus the function Q is a CBF, hence 
by Theorem 12.81 i^ipY Ip = pp/Q{p) is a CBF. □ 

In the case of power law attenuation = C p" the function k{pY/p = 
B^p+2BCp''+C^p'^"-^ isaCBFif andonlyif 1/2 < a < 1. If a < l/2then 
the exponent 2a — 1 is negative and the function Pip^/p is not monotone, 
hence it is not a BF; consequently the relaxation modulus is not completely 
monotonic. The function f{p) := p/{l+p) is not a possible candidate for the 
wave number function because p^ / f{pY is not a Bernstein function. Other 
counterexamples of this kind can be found in jof. 

We now examine the class of functions K{p) such that i^{pY/p G 0. 
Eq. ( ITOl) implies that k{pY /p G if and only if G and (3{p)'^/p G 0. 
We shall use the following notation: 

0" := {r I / G 0}, « G M 

For any real function / defined on R+ 

/0:= {fg\geil} 

where fg denotes the pointwise product of functions. 

Theorem 6.2. Q{p) G z/ and only z/ k G fl p^^"^ 0. 

Proof. If the function Q{p) = p^ / K^pY is a CBF then the function Q{pY^'^ is 
also a CBF and the complex wave number function K{p) = p^^'^ p/lQip)]^^"^ 
is therefore also a CBF. On the other hand k{p) = p^/^p^/^ [p/Qip)Y^'^- The 
third factor is the square root of a CBF and hence a CBF itself. Hence 
K, G p^/^0 are also CBFs. We have thus proved the "only if part. 

For the converse we shall use the following identity (Theorem 7.11 in (lol |) 
for a G [-1, 1]: 

0" = 0np"-^0 (15) 

i.e. 

0" = {/G0b^-"/(p)G0} 
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Hence p^/^ H^/^ = p^/^O fl O.. By our hypothesis k E p^/^H^/^, therefore 
= p f{p), where f E Q.. Thus Q{p) = p"^ / K{pY = p/ fijp) E O.. □ 

We have thus proved that the mapping Q — )■ k is a bijective mapping of 
onto the space K, := Hflp^/^ 0. This fact imphes that the wave number func- 
tion has another integral representation, viz. k(p) = p^^"^ J^^ oo[(^'~'~'")~^ A(dr) 
with A e 971. 



7. Finite propagation speed. 

Since 



son^ \ f \P\ +r^P 
^Q{p)= — — —^fi{dr 

i[o,oo[ \p + r\^ 



'[0,oo[ 

the function \Q{p) \ is a non- decreasing function of |p| in the right half complex 
p-plane. It increases to infinity for \p\ — )■ oo in the right half complex p-plane. 
If the spectral measure u has infinite mass then the real part of the function 
K,{p) = p/coo + P{p) increases to infinity in the right half complex p-plane 
(Theorem 15.31) . 



Lemma 7.1. Let Coo < oo. 

The function f{p) := [Q{p) k{p)]^^ e^^^'^^^ is analytic in the right half 
plane and tends to for \p\ — oo uniformly with respect to argp E [—n/2, 7i/2] . 

Proof By Theorem O | / (p) | < \l/[Q{p) n{p)]\ in C+. If 3fJp > then 



in view of the inequality |p + r | < r + Let a > be sufficiently large so 
that /i([0,a]) > 0. The integrand of the last integral is non-negative, hence 

For an arbitrary positive Ei < /i([0,a]) a sufficiently large Ri can be found 
so that ^Q{p) > /i([0,a]) — ei and \Q{p)\ > /i([0,a]) — ei for \p\ > Ri. 

In view of the inequality = l/{\p\ \B + /3{p)/p\) Theorem 15.41 

implies that for every positive e < 1 there is a positive R such that < 
[1/(1 - e)]/\p\ for all p with \p\ > R and -n/2 < argp < n/2. □ 
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Theorem 7.2. // Cqo < oo, then the Green's function vanishes for t < 

/ ^oo • 

Proof. Let t < x/Cqo- 

The Green's function is given by an expression of the form u^^\t.,x) = 
^ Jb /(p)e^'^*~^/'^°°''dp, where B is the Bromwich contour running from 77 — ioo 
to ?7 + ioo, 7] > 0, parallel to the imaginary axis, A is a constant and the 
function / is defined in Lemma I7.1[ Consider the complex contour Cr : p = 
re^'^, with running from 7i/2 to —n/2 and a fixed r > 0. Lemma [7.11 and 
Jordan's lemma imply that the integral A f{p) exp(p [t — x/coo))dp tends 
to as r —7- 0. Let Br be the straight line contour running from 77 — ir to 
r] + ir. By the Cauchy theorem the integral of f{p) eP^^~^/^°°) over the contour 
Br+Cr vanishes. Hence, taking the limit r — )■ 00, the integral over B vanishes. 

Consequently u^^\t,x) = for t < x/coo- □ 

Corollary 7.3. If c^o < 00 then u^^\t,x) vanishes for t < |x|/coo- 

Proof. The thesis follows from eq. ([7]). □ 



8. Dependence of wavefront smoothing on the spectral density. 

A frequent feature of wave propagation in real viscoelastic media is wave- 
front smoothing 18l-l22j. Using a terminology often adopted in mechanics, 
many linear viscoelastic media do not allow non-trivial discontinuity waves. 
It will now be shown that absence or presence of non-trivial discontinuity 
waves depends on a property of the dispersion-attenuation spectral measure. 
Wavefront smoothing is ultimately due to the singularity of the derivative 
of the relaxation modulus at 0. We shall now relate it to the asymptotic 
properties of the attenuation-dispersion spectral measure. 

If all the moments of the spectral density function h 

h{r) dr 

Jo 

are finite, then 



00 



m=P -oo }_^(-l)"a„p-" (16) 

Jo p + r ^ 



n=0 
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16l . Il7j . The dominating term is a positive constant Oq. A special case is a 
finite bandwidth spectral density such as h{r) = Kx[a,b]{r), K > 0, < a < 
b<oo. In this case (3{p) = K \n[{p + b)/{p + a)] r^oo K {b - a) + O [p'^]. 

If the function h decays at an algebraic rate then some higher order 
moments are infinite and f|T6l) does not hold. Assuming the asymptotic ex- 
pansion of the spectral density 

oo 

^(r)~oo5^&nr-"-", 0<a<l (17) 

n=0 

the dissipation-attenuation function has the following asymptotic expansion 



at infinity [17 



N-l 



m ~oo - nc„p-"] + Rm (18) 

where 



n=0 



/>oo " 

/ /„(r)r"-Mt, f^{r):=h{r)-y^b 
Jo 



k=0 

P Jo P + r 

The asymptotic expansion of /3 is valid in the entire cut complex plane, 
I argpl < 71 [L7\. The first term of the expansion of the dissipation-attenuation 
function (3 is now bo while the attenuation function A{u) := 3?/3(— icu) ~oo 
bo sin(7ra;) u^'". 

Consequently, if the spectral density h decays algebraically at infinity, 
then I exp(— ia;(t — x/coo) — /3(— iw) |a;|)| < | exp(— ^(cj) |a;|)| vanishes asymp- 
totically like exp (— 6o sin(7ra) o;^"" r). Hence the integral ([8]) is absolutely 
convergent and therefore the function u^^^ is continuous. Furthermore, the 
derivatives of D'^D'^u^^^t, x) are also given by absolutely convergent inte- 
grals. Consequently u^^\t,x) is a smooth function of both arguments in 
]R_i_ X ]R_|_. In particular it has continuous derivatives of arbitrary order at the 
wavefront t = x/coo- Thus the wavefront does not carry any discontinuity of 
the Green's function nor its derivatives. Since the Green's function vanishes 
for t < a;/coo, it gradually decays to with all its derivatives. 

Eq. ([7]) implies that u^^\t,x) is also a smooth function of (t,x) G M+ x 
(M3\{0}). 
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The case of a strongly singular relaxation modulus is analyzed in Theo- 
rem [STTl 



Theorem 8.1. If G{t) = /j^^je-'Xdr) with e Tl and At([0,r]) 
r"l{r), where < a < 1 and I is slowly varying at infinity, then Gq = oo, 

Coo = OO, Q{p) C,P"/(P), <P) ~oo p^/2c-V2p7/(p)-l/2^ 1/2 < 7 < 1 

and udO^r]) ~oo p^^'^ Ca^^'^ l{p)~^^'^ , whereby := 1— a/2, := Tta/ sin(7ra). 

If /i([0,r]) ~oo r°'l{r), < a < 1, then, by the Karamata Abelian The- 
orem (Theorem I B.Sp . G{t) ~o t~"/(l/t)/r(l — a). In this case the Green's 
function is a complex analytic function of {t, x) in a neighborhood of M+ x M 



18|. (This follows from a more general fact that the Green's function is 
analytic except at the wavefront; in the strongly singular case the Green's 
function does not have a wavefront). Analyticity implies that the Green's 
function cannot vanish on any open subset of M+ x M unless it is identically 
zero. 

The case of weakly singular relaxation modulus is somewhat more com- 
plicated. The relaxation modulus is weakly singular if G is bounded and 

G'{t) bt-"l{t), 6 > 0, < a < 1. In this case /(r) := p(\0,r[) ~oo 

//(]0,oo[) — ar~", with a = bT{l — a), Gq = /i({0}) < oo, Go — Goo = 
yu(]0,oo[). / is a function of slow variation at infinity vanishing at 0. By 
Valiron's Theorem 

Q{p) = G^+p [ ~oo Goo + /i(]0, oo[) - ap-" = Go - ap-" 

i]o,oo[ P + r 

Hence 

/ \ 1/2 

l^ip) ~oo ( 7^ ) ^ ,^ ,1/2 

This implies that /3(p) ~oo [ct/(2Go Coo)] p'^ where 7 := 1 — a G ]0, 1[, and 

/r IN a sin(7ra) , „ 

2Go Coo 7r(l - a) 

The dissipation-attenuation spectral density has an infinite bandwidth and 
an algebraic decay at infinity and the Green's functions are infinitely smooth 
at the wavefront x = Coot. 
Summarizing, 
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(i) if z/(dr) = h{r) dr and all the moments of the spectral density h are 
finite, then the Green's function can have discontinuities at the wave- 
front; 

(i) if z/(dr) = h{r) dr and the spectral density function h decays at an 
algebraic rate then the wavefront does not carry any discontinuity of 
the Green's function nor its derivatives of arbitrary order. 



9. Attenuation and dispersion functions. 

Define the attenuation function A and the dispersion function V by the 
equations 

A{u) := ^f3{-iu) = ^K{-iuj) [ 4^ (19) 

i]o,oo[ + 

V{u) := -$>/3(-ic.) =uj '-^^ (20) 
By Theorem 15.41 

^(tj) = Ooo[w], V{uj)=o^[uj] (21) 

Note that Aiu) > and sgn Vlu) = sgn u. Recalling the Green's function 
(jH]) and the definition of k(p), this implies an outgoing sense of propagation 
and a non-negative attenuation along each radial direction. 

Since x/{l + x) is an increasing function, eq. ( |T9l) implies that the atten- 
uation function A is non-decreasing. 

The phase speed c{u) := —uf^n^—iuj) = l/(l/coo — 3'/3(— iw)/^) is re- 
lated to the dispersion function by the equation 

+ -^-^ (22) 



C{u) Coo W 

Note that c{co) < Coo because V^u) > 0. Moreover V^u) /i^ = ijo oo[ z/(dr) / (c<;^+ 
r^) is a non-increasing function of u. Therefore the phase speed c{u) is a 
non-decreasing function of frequency. 
Eq. (I22]) and eq. (El]) imply that 

lim c(u;) = Coo (23) 
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If Q{0) = Goo > then, in view of eq. (Q, \imp_:,Q k{p) / p = {p/GooY^'^ 
for all p G C. The inverse of this limit will be denoted by the symbol Cq. In 
particular, for p = —iu, w G M, this proves the following theorem: 

Theorem 9.1. //Goo > then 

lim ^(a;)/a; = 0; D := \imV(u)/u < oo 

and 1/co = lim^_^o[l/c(w)] = l/coo + D. 

Note that D > and D > unless G{t) = const. Under the hypotheses 
of Theorem 19. II the phase speed satisfies the inequality Cq < c{u) < Cqo, while 
< c(w) < Coo in the opposite case. 

Corollary 9.2. Goo > entails that 

D=[ ■^<o. (24) 

Proof. Indeed, 

V{uj) f z/(dr) 1 



For w — )■ the function l/(l+ci;^/r^) increases monotonically to 1. Therefore, 
if it is assumed that the integral in fl2^ is infinite then, by the Fatou lemma, 
V{uj)/uj tends to infinity for w — ?■ 0. On the other hand, if the inequality in 
is satisfied, then 

Z,= ,ta5M=/ !1M (26) 

by the Lebesgue Dominated Convergence Theorem. □ 

10. High-frequency behavior and the Kramers-Kronig dispersion 
relations. 

The Kramers-Kronig (K-K) dispersion relations are the Sochocki-Plemelj 
formulae following from the fact the dissipation- attenuation function P{p) is 
the Laplace transform of a causal distribution F{t). Indeed, 
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where L denotes the Laplace transformation. The integral representation 
of the function /3 shows that all the singularities of the complex analytic 
function P{p)/p lie in the closure of the left half p-plane. By Jordan's lemma 
and Theorem 15.41 the source function g of Pi{p) := P{p)/p vanishes for t < 0. 
Applying the inverse Laplace transformation L"^ to both sides of eq. (|26|1 we 
have 

g{t) = L-\f3,)it) = L{u) = [ e-^'u{dr) 

J[0,oo[ 

for t > with u G VJt. By Theorem 12.41 the function g is LICM. Hence its 
primitive 



fit) :-- 



[ g{s)ds= [ [l-e~*"] z/(dr 

Jo J]0,oo[ 



is a Bernstein function continuous over [0, oo[ and /(O) = 0. It follows that 
the dissipation- attenuation function P{p) is the Laplace transform of the 
causal distribution f of second order. For f3{p) = C p"', with < a < 1, 
C > 0, the function / is Ctl-''/T{2 - a). 

Theorem 10.1. If the v measure oj the closed segment [0,r] is a function 
of r regularly varying with index 1 — a, 

z/([0,r])~oo r^-°/(r) 
where I is a function of slow variation at infinity, then a > and 

A{u)^^,p^^^y--l{u) (27) 
2 cos{aTT/2j 

V{uj)^^SLz^^uj'-l{uj) (28) 
2 sm(a7r/z) 

Proof. Since u G eq. ([2]) implies that a > 0. 

LetF(u;) := j^^^^^{uj^+r^)-^ p{dr) = /^^^^(w^^r^)-! dz/([0, r]), where the 
integral on the extreme right-hand side is a Stieltjes integral. The function F 
can be expressed in terms of a Stieltjes transform by changing the integration 
variable. In terms of a new measure /x([0, s\) := z/([0, y/s\) = s(i~")/2 / (y/s), 

[0,oo[ ■ 



FM=/ ^ (29) 
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By Valiron's theorem (Theorem I B.4p 

^ ' 2cos(7ra/2) ^ ' 

This proves the theorem. □ 

The function F in eq. (129|) is differentiable for a; > 0, hence the atten- 
uation and dispersion functions A and V are differentiable. The function 
00^'^ A{uj) ~oo by Theorem 110. ![ hence it belongs to the space V'^^^ 
(eq. (1.8.12) in [23| Sec. 1.7-1.8). Consequently under the above hypothe- 
ses the functions A{ijj) = 3ft/3(— iw) and T>{uj) = — 3/3(— iw) satisfy the K-K 
dispersion relations with one subtraction: 

^(,) _ AM = -(^vp r (30) 



7^ J-oc i^' - ^o) (W - W) 

where "vp" indicates that the integral is to be taken in the sense of principal 
value. 



11. Attenuation and dispersion functions. Low- frequency behav- 
ior. 

Low-frequency behavior of the attenuation function provides a useful test 
whether a specimen of the material subjected to constant strain for t > 
relaxes to zero stress - that is, whether Goo = 0. Materials with vanishing 
Goo are known as viscoelastic fluids. Scalar models considered in this paper 
represent either longitudinal or shear waves. It turns out that the same 
material can behave under tension or compression like a viscoelastic solid 
and under shear strain as a viscoelastic fluid. 

If z/ G 971 and the u measure z/([0,r]) of the closed segment [0,r] is a 
function of r regularly varying at with index 7, then 7 > 0. 

Theorem 11.1. // the Radon measure u & DJl is regularly varying at with 
i/([0,r]) = l{r), where /(r) is slowly varying at and < 7 < 1, then 

(i) A is regularly varying at with 

2 sm(77t/2) 

with I slowly varying at 0; 
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(ii) V is regularly varying at with 



cos(77t/2) 

(iii) the Q factor is asymptotically constant at 0, Q{uj) ~o 47t cot (771/2). 
Proof. Ad (i) 



where the integration variable has been changed to s = and /i([0,s]) ~o 
s'^/^/(l/s). For example, if z/(dr) = h(r) dr with h{r) ~o 'jr"'~^l(r), then 
Ai(ds) = g{s)ds with g{s) = h{^/s) /{2^) ~o {l/2)s^/^-H{^). The func- 
tion l{^/s) is slowly varying at 0, hence Lemma [ B.4I implies the thesis. 
Ad (ii) 

^/ N f ruidr) f ^/s n(ds) 



Valiron's Theorem implies the thesis. □ 

Theorem II 1.1 1 is applicable to viscoelastic fluids. In the case of viscoelastic 
solids eq. flMj) implies that 7 > 1 and the hypotheses of Theorem 111.11 are not 
satisfied. In this case Theorem 19.11 implies that Viuj) = D oj + Oq[u]. If A{uj) 
is regularly varying at then A{(jj) ~o \^\°'^^ K'^) ■> where / is slowly varying 
at and either a > or /(O) = 0, for example A{(jj) ~o C \uj\ \'n.^{l + bj) with 
C,/3 > 0. 

Attenuation of longitudinal waves in polymers and bio-tissues exhibits the 
power law behavior A{uj) ~o C|a;|^"'"°'. The asymptotic behavior described 
by (i) of Theorem 111.11 is observed for shear wave attenuation in minerals 
(e.g. H). 



12. A few examples of analytic viscoelastic models. 

We shall now demonstrate an application of the spectral method to the 
numerical determination of attenuation and dispersion in viscoelastic media 
defined by the most popular analytic expressions for the complex modulus. 
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12.1. Power law. 

If z/([0,r]) = ar''', < 7 < 1, then A{u}) = Au'^, where 



A = aj 



12.2. Finite bandwidth. 

If i/(dr) = C X[a,b]ij') dr, with 0<a<6<oo, C>0, then the attenuation 
function 

A{uj) = Ccu [tan~^(6/a;) - tan"^(a/a;)] 
is asymptotically constant at infinity and O [o;^] at 0: 

'oo C [b- a) 

On the other hand I/cq = limp^QK{p)/p = 1/coo + C ln(6/a), hence T>{u) ~o 
C \n{b/a) u. 

12.3. The Cole- Cole relaxation model. 

The Cole-Cole model is defined by the Laplace transform of the relaxation 
modulus of the form 

gg"(p) = ^ W(S)r ^ a>l, 0<a<l, G^>0, r>0 

„ (32) 



25|, |26[. The Cole-Cole model [29[ and its generalizations [27[ are equivalent 



to linear viscoelastic models based on fractional derivatives considered in |28| . 
The function G^*" is bounded completely monotonic (Appendix 1"^ . Using 
the formula 

^*e-^*E„,,(-Ar)dt = ^^ (33) 

30l | the relaxation modulus can be expressed in terms of the Mittag- 
LefHer function Ea = 

G'S''(t) = Goo 9{t) [1 + (a - 1) i-r/b)] (34) 

The relaxation modulus G^*" has a finite limit Go = Goo a at t = 0. 
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The complex wave number function of the Cole-Cole model is given by 
the formula 

1 + {rpr ^ 



Co 



1 + a{TpY 



Bp + p 



z/(dr) 

0,oo[ P + 



(35) 



An explicit integral expression for the attenuation and dispersion functions 
of the Cole-Cole model can be derived from eqs (fTO] - [TTl) and (fT^ . In the first 
place Corollary 15.51 implies that 



B 



lim 



-1/2 



I Co 



(36) 



p 

The wavefront speed equals Coo = 1/B. 

The spectral density /i^*" of the Cole-Cole model is given by the formula 



Let 



-53 



71 Co 



(r exp(i7r)) 
r exp (in) 

1 + {tpY 



1/2 



n := SR- 



1 r 1 + (rr)" exp(ia7r) 

^ 

nco 1 + a (rr)° exp(iQ;7r) 
1 + (rp)" 



1/2 



1-Fa(rp)°' ■ l + a(rp)" 

where p = r exp(i7t). The imaginary part Y of [(1 + (tjo)") / (1 + a (rp)")]"*^^^ 
is a solution of the bi-quadratic equation + TZY^ — J'^/A = and is 

non-negative, hence Y = ^/VWTj^^^/V^. Now J = Ji/Z^, 7^ = 



TZi/Z'^, where Ji 



1) sin(7ra) (rr)°, 7^l 



1 + g ( rr)^" + [a + 

1) cos(7ra;) (rr)" and Z = a/1 + (rr)^"' + 2acos(7rQ;) (rr)". Hence the 
dispersion-attenuation spectral measure z/ of the Cole-Cole model has a den- 
sity 



h 



CC/ 



7^l| 



(Ji/7^l)2 



sgn(7^i 



71 Co 



Note that 



h 



CC/ 



fei/^ g-l 

TtCoo a 

1 /'^ 



sm(Q;7tj (rrj 
- 1) sin(a;7r) (rr) 



(37) 



(38) 



and the function h'^^{r)/r is integrable. Thus T>^'~^{ijj) ~o Z^w. The asymp- 
totic properties of the attenuation function follow from Theorem 111.11 



iSsin(W2) (™ 



,1-a 



CO 



- 1) sin(a7r/2) (™)i+" 
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(39) 



For a — )• 1 the dispersion-attenuation spectral density h^^ tends to the 
spectral density of the Standard Linear Solid 



^ I 0, 7^1 > 

V I ^ < ^ ^ 

Ttco \a r—l\ ' 



or, equivalently 



,SLS' > I ± TV 



r 



TtCo V O'T"?" 



^X[i/a,i](rr) (41) 



In contrast to the Cole-Cole model, the spectrum of the SLS has a finite 
bandwidth and h^^^{r) dr < oo. 

The attenuation functions of the Cole-Cole and the SLS are now given 
by the explicit expressions 

Jo + 

for < q; < 1 and a = 1, respectively. The dispersion functions are given by 

< a < 1. The attenuation function and the phase speed c{u) with Cqo = 
5000 m/s are shown in Fig. [1] 

Equations ( l42l) and ( l43l) are convenient for numerical computation of the 
attenuation and dispersion functions. 

12.4- The Havriliak-Negami relaxation. 

The Havriliak-Negami relaxation model is defined by the Laplace trans- 
form of the relaxation modulus of the form 



GS(P) = ^ 

' ' p 



1 



0<6<1, 0<a<l, 0<7<1, Go>0 (44) 



3ll |. The corresponding relaxation modulus is given by a 4-parameter for- 
mula involving the Prabhakar Mittag-LefHer function ^-y 0] • The Havriliak- 
Negami relaxation modulus G™ is a CM function (Appendix 1"^ . 
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(a) Attenuation. (b) Phase speed. 

Figure 1 : Attenuation function and phase speed plotted vs logarithm of angular frequency 
in MHz for the SLS and Cole-Cole models, a = 1.5, Coo = 5km/s, t — 10~^^s. 
Solid line: Cole-Cole, a = 1/2; dashed line: Standard Linear Solid. 



For 7 = 1 the Havriliak-Negami relaxation modulus reduces to the Cole- 
Cole model with a = 1/(1 — 6). 

Eq. © implies that c^o = {G^/ pf^ and = [1 - 6/ (1 -f {rpYfY^'^- 

1/coo- Hence the spectral density is given by the expression 

TTCoo 

where Z := \ — h/Y and Y := {1 + (rr)" exp(i7rQ; ))'''. 
Substitution of the expressions 

|r| = g(r) ■= (1 + 2(rr)" cos(7ra) + (rrfY^^ 

3ftF = g{r) cos(7/(r)), QY = g{r) sin(7/(r)) 

where 

/(r) := tan""*^ ((rr)" sin(7r«)/ (1 -|- (rr)" cos(7ra))) 
in the formulae 

spy cky 

yields the formula 

\Z\ = k{r) := ^Jl + / girf - 2b cos(7/(r))/^(r) 
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Hence, using the identity 



^zV2 J\z\-m 



and noting that the lower sign apphes in the problem at hand, the final result 



IS 



/,HN/ N _ 1 v/A;(r) - l + &cos(7/(r))/^(r) 

Note that /(O) = and therefore fc(0) = 1 — 6/(7(0). This implies that 
/,HN(0) = and D = /o"[/i^,^(r)/r]dr < 00. 

In comparison with the Cole-Cole model the parameter 7 in the Havriliak- 
Negami model allows decoupling the high-frequency asymptotics of the wave 
number function from its low-frequency asymptotics. The high-frequency 
asymptotics Q{j)) ~oo (1 — 6 (tp)""'^) implies that = p/coo+/3(p) ~oo 
p/coo + p [1 + b{Tp)~"'^] and thus f3{p) ~oo (6/coo T)(rp)^~"'^. On the other 
hand Q{p) ~o Goo [1 + 67/ (1 — 6 (tjo)"] and therefore k{p) ~o p/cq x 
{1 — 67 (rp)°"'"^/[2(l — 6) r]}. The lowest frequencies propagate with speeds 
close to Co and exhibit the attenuation A^u) ~o ^7 sin(a7t/2)(ru;)^+"/[2(l — 
6)rco], in agreement with experiments. 

Due to its fiexibility achieved by a minimal number of parameters and 
asymmetric shape of the peak of the loss modulus the Havriliak-Negami 
relaxation model is frequently used in modeling the alpha relaxation in the 



mechanical and dielectric data in polymers and glass-forming liquids |32H35 



12.5. The Cole-Davidson relaxation. 

The Cole-Davidson relaxation modulus is given by (jHj) with a = 1. The 
spectral density h^^ of the Cole-Davidson model cannot however be obtained 
by substituting a = 1 in eq. f H5|) . For a = 1 the branch cut in the complex 
p-plane is reduced to the segment ] — 00, — 1/r] and 



Y 



(rr — exp(i7r7), rr > 1 
(r r — 1)'^, r r < 1 



Hence 



, 6 cos(7r7) (rr — 1) rr > 1 ) 6 sin(7t7) (rr — 1) rr > 1 

utZ = 1 — < , isZ — 



6(1 — rr) rr < 1 10, rr < 1 
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(a) Attenuation function. (b) Phase speed. 

Figure 2: Plots of attenuation and phase speed vs logarithm of angular frequency in MHz 
for the Cole-Cole, Havriliak-Negami and Cole-Davidson models, b — 0.5, Coo — 5km/s, 
T = 10"" s. 

Solid line: Cole-Davidson, /3 = 1/2; dashed line: Havriliak-Negami, a = 1/1.3, (3 = 1-3/2; 
dot-dashed: Havriliak-Negami, a = 1.3/2, /3 = 1/1.3; dotted: Cole-Cole, a = 1/2. 



Consequently 



/i?°(r) = ^ ^'^^^'^ / ' "^""^'^'^ ^(rr- 1) (46) 



1 ^/ ki{r) — 1 + 6003(717) 
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with fci(r) := a/1 + b^/\l - rr |2t - 26 003(717)/! 1 - Tr\^. 

The low-frequency asymptotics of the Cole-Davidson attenuation function 
can be calculated by substituting eq. (jH]) in Qj and using the definition ([T9|) 
of A yields a behavior characteristic of a spectrum which does not extend to 

0: A{u) ~o7&^c^V[2co(l-&)] (dSec.MM- 

Ultrasonic data for polymers and bio-tissues cover the range 0-250 MHz, 
which is much below the inverse of the characteristic relaxation time 1/r = 
10^ MHz for polymers. Therefore they are expected to be matched by the 
low-frequency asymptotics of the attenuation function and phase speed. At- 
tenuation and phase speed for the Cole-Cole, Havriliak-Negami and Cole- 
Davidson models in the low-frequency range are shown in Fig. [2J The plots 
have been obtained by numerical integration of eqs f|T9] - [20|) using eqs psj) . 
fH6|) and fl22l) . The asymptotic high-frequency exponent equals 1/2 for all 
these cases. 
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13. Plane waves and minimum phase signals. 

The K-K dispersion relations are closely connected to the minimum phase 
properties of the propagators. 

Definition 13.1. A causal tempered distribution T is said to be minimum 
phase if its Fourier transform T and its inverse 1/T{uj) are analytic in the 
upper half-plane C+. 

The first condition implies that the poles and branch cuts of T lie in the 
lower half u plane. The second condition implies that T does not vanish in 
the upper half plane C"*". 

A shifted function or distribution Trf is defined by the formula {Trf){t) = 
f{t — t). Its Fourier transform 7^/ (a;) = /(w) e''^'^ has the same zeros and 
singularities as f{uj). 

The following proposition is obvious. 

Proposition 13.2. 

(i) // the tempered distributions T and S are minimum phase then their 
convolution is a minimum phase distribution. 

(ii) // the convolution of the tempered distribution S with an arbitrary min- 
imum phase tempered distribution is a minimum phase tempered distri- 
bution, then S is a minimum phase tempered distribution. 

(iii) // the distribution T is minimum phase, then the shifted distribution 
TrT is minimum phase distribution. 

The convolution T*S can be viewed as an application of a linear operator 
T to a signal S. A minimum phase operator is defined as a linear operator 
that preserves the class of minimum phase tempered distributions. It follows 
from Proposition 113.21 that a minimum phase tempered distribution M can 
be viewed either as a minimum phase signal or a minimum phase linear 
convolution operator (filter). 

Theorem 13.3. A tempered distribution T is minimum phase if and only if 
ln(T) is analytic in the upper half plane C^. 

Proof. ln(T) is analytic in C+ if and only if T is analytic in and does 
not vanish there. On the other hand, if T is analytic in C^, then 1/T is 
analytic in C+ if and only if T does not vanish there. The thesis follows by 
comparison of the last two statements. □ 
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Let the distribution or function / be causal. Its Fourier transform / is an 
analytic function in C^. Assume that ^f{u) — oo for — )■ oo in w G C"*" 
is an increasing function of |a;| in C"*", with 3fJ/(a;) = O [\u}\'^] for some 7 > 0. 
The functions e^-'^^'^-* are analytic in and the function e"'^^'^-* tends to zero 
as — > 00 in C_|_. The inverse Fourier transform 
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00 



—00 



is absolutely convergent and hence a continuous function of t. By Jordan's 
lemma the function g{t) is causal. It is also a minimum phase distribution. 

Theorem 13.4. Let T be a minimum phase tempered distribution. Then 
ln(T) satisfies the K-K dispersion relations with two subtractions. 

Proof, (sketchy) By Theorem 113.31 ln(T) is an analytic function in C^. Since 
T is causal, by the Paley- Wiener Theorem (theorem I A.ip \n{T{uj)) = Ooo[i^]- 
The function H{uj) := \n{T{u)))/uj'^ is analytic in C"*" and it uniformly tends 
to as a; tends to 00. Hence H is the Fourier transform of a causal function 
h{t). Since H{u) = o the inverse Fourier transform of H is absolutely 

convergent, the function f & C^. Hence / is a distribution of second order 
and ln(T) satisfies the K-K dispersion relations with two subtractions. □ 

The Fourier transform of the Green's function u^^\t,x) is the product of 
three factors 

^^^^ ^ 2k{' iu) e-^(-^^)^ (47) 

The inverse Fourier transform of fBTl) is a convolution 

u{t, x) := 5{t - x/coo) *t M{t, x) *t s{t) = M{t - x/c^,x) n Sit) (48) 

where S{t) is the inverse Fourier transform of the first factor and M{-,x) is 
a minimum phase signal for each x > 0. The two first factors can be viewed 
as transfer functions or filters applied to the source signal. 

Since K,{p) is a CBF function, the singularities of the function k{—iuj) = 
iK(— ia;) lie on the negative imaginary axis. By Corollary I C. 51 its zeros also 
lie on the negative imaginary axis. Hence the first factor is a minimum phase 
filter. 

The third factor has two crucial properties 
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1. it is the Fourier transform of a causal tempered distribution M; 

2. the real and imaginary parts of the phase ~f3{—iu)x of the Fourier 
transform M of M satisfy the K-K dispersion relations. 

These two properties identify M[t) as a minimum phase signal. Note that by 
definition a minimum phase signal has an onset at zero time. In seismology 
and acoustics the onset time of the minimum phase signal can be arbitrary. 
The signal in (jlSj) is shifted and it starts at t = x/coo, the travel time to the 
point X. 

The second factor in ( H7|) is an all-pass filter because it has a constant 
unit amplitude. It is a shift operator acting on the third factor. A shifted 
minimum phase function is again a minimum phase function. 

Hence u^^\-,x) is a minimum phase signal. If the source function s{t) is 
a minimum phase signal, then the convolution s *t u^'^\-,x) is a minimum 
phase signal for every x. 



14. Concluding remarks. 

We have demonstrated the utility of the spectral decomposition of the 
wave number function for investigating general properties of the attenuation 
and dispersion in a viscoelastic medium with a positive relaxation spectrum. 
It has also been shown that spectral decomposition of attenuation provides 
a very efficient method for numerical computation of attenuation in some 
popular viscoelastic models of polymers and soft matter. 

Theorem 15.11 implies that the complex analytic continuation A{u) of the 
attenuation function, defined on the complex plane except on the negative 



imaginary axis, is a Herglotz function. Weaver and Pao [36| derived the K-K 
dispersion relations from this property. The derivation of the K-K dispersion 
relations adopted in this paper is based on a constitutive assumption. Weaver 
and Pao's assumption is not entirely based on a physical argument because 
it concerns the wave number function at complex-valued frequencies. 
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A. The Paley- Wiener Theorem 

How can one verify whether a candidate function M{uj) for the complex 
modulus is the Fourier transform of a function vanishing for t < 0? 

The best available answer to this question is Theorem XII in Paley and 
Wiener (37] : 

Theorem A.l. Let the function M{u) be square integrable. 
If 

|ln|M(w)|| 



1 + ^2 



■ dcj < oo, 



then there exists a function G : M — )■ C with support in the closed positive 
real semi-axis M+ such that M is the Fourier transform of G. The converse 
of this implication is also true. 

Note that the hypotheses of the theorem involve conditions on the abso- 
lute value of M{u) only. 



B. Regular variation. Valiron's and Karamata's theorems. 

Definition B.l. A real function f is said to be regularly varying at a, where 
a = or oo, if limx^ a f i^x) / f (x) is finite. 

f is said to be rapidly varying with index oo if 



lim f{Xx)/f{x) 
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A real function f on is said to be rapidly varying at a, where a = 
or oo, with index — oo if 



Theorem B.2. // / is regularly varying at a, where a = or oo, then 
lima,_j.(j /(Ax)//(x) = A", where a is a real number. 

The number a is called the index of the function / at a. 

Definition B.3. A real function f on M+ is said to be slowly varying at a 
if it is regularly varying at a with index 0. 

A regularly varying function / with index a at a can be expressed as the 
product l{x), where / is slowly varying at a. 

Examples: ln(H-a;) is slowly varying at infinity; + is slowly varying 
at if a < and regularly varying at infinity with index a if a > 0; exp(— x) 
is rapidly varying at infinity with index — oo. 

Theorem B.4. (Valiron 1911, see iSdJ ) If F is an increasing function 
satisfying the condition limt^o- -^(^) = ^'^^ ^he function (j) is given by the 
Stieltjes integral 



where < A < 1 and the function I is slowly varying at infinity. 

Valiron 's Theorem also holds when oo is replaced by 0. 

Theorem B.5. (Feller fMj, Chapter XIII; Seneta fJ^J, Theorem 2.3) 
If F is non- decreasing right- continuous, I is slowly varying at 0, 'y > 0, and 
the Laplace- Stieltjes transform (f){p) of F converges for p > a, where a is a 
positive real number, then 




A > 1 



A < 1 




oo 



r/(iA)/r(i + 7) 
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C. Stieltjes functions. Zeros of a CBF. 

Definition C.l. A n on- negative real function f on with the integral 
representation 

b f fi{dr) 



f{x) = a+-+ ^ (C.l) 

X i]o,oo[ x + r 

where a,b > and fi G QJt, is called a Stieltjes function. 

Eq. (IC.ip can also be expressed in the form /(x) = a + f^Q^^fi{dr)/{x + r) 
with Ai({0}) = b. 

A Stieltjes function f{x) has a complex analytic continuation f{z) in the 
complex plane cut along the negative real axis, | arg(2;) | < n. The function 
f{z) has the integral representation ( IC.lll with x replaced by z. It follows 
immediately that f{z) maps to and vice versa. 

Lemma C.2. A Stieltjes function has a finite value at every point of the 
complex plane outside the closure of the negative real axis. 

Proof. Let z = x + iy. 

The inequality I2: + t| = A/(r + x)^ + y'^ > r + x implies for a; > the 
inequality 



/i(dr) 

[o,oo[ z + r 



< I ^<oo 

[o,oo[ x + r 



(cf Remark [T] ). 

U y > then we can find such a number -i? that < < 1 and q :- 
"i^x + (1 — {})y > 0. Hence 



V(r + x)2 + ?/2 y ^(^r + x) + (l - 1!}) y > 1:} {r + qf-d) 

Hence 



/i(dr) , 1 /■ /x(dr) 

< -r / r-r < OO 



[0,00 



□ 



Theorem C.3. ^ 

yln analytic function f regular in C+ and satisfying the inequality (i) Qf{z) < 
in C"*" and (ii) f{x) > for x G IR+ can be expressed in the form 



f\z) = fo + 
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[o,oo[ y + z 



where /o = lim2_j.oo f{z) and fi G OJt. The limit lim2_j.oo f{z) exists and equals 
/q. The measure fi is uniquely defined as the weak limit 

0(y)/i(dy) = -- lim / ^(y)Qf(-y + ie)dy (C.2) 

[0,oo[ 71 e^O+J[o,oo[ 

for every continuous function (p with compact support. 

Theorem C.4. The inverse 1/f of a non-zero CBF f is a Stieltjes function. 
Proof. The integral representation (j4]) imphes that 

f{l/z) = ai + b-+ / = ai + b- 



z 



0,oo[ 

^ + zy z J]o,oo[ u + y 



where oi := a + p({0}) > 0, m = 1/y. Hence the function f{l/z) satisfies (i) 
and (ii) of Theorem I C.3I In particular, (ii) implies that '^f{l/z)~^ > for 
z e C"*". Substitute z = Clearly, sgnS^C = — sgn^Sz and, by reflection, 
^fiCy^ < for C G C+. Since /(x)"^ > as weU, /(x)"^ is a Stieltjes 
function. □ 

Corollary C.5. A CBF does not vanish outside the closed negative real 
half- axis. 

Proof. Let / be a non-zero CBF. If /(a) = then the Stieltjes function 1/f 
has a pole at a. Hence a G M_. □ 

D. Proof of Theorem Em 

Proof. By the Fubini theorem 

1 r r ri r i _ q-t 



i[0,oo[ J[0,oo[ Jo i[0,oo[ ^ 

If / if integrable over [0,1] then the last integral is convergent. Since the 
integrand is non-negative and the Radon measure n is positive, 

1 - e"*" 

/i(dr) < oo 

l,oo[ 



But 

J[l,oo[ ^ ^[l,oo[ ^ J[l,oo[^ + ^ 
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hence the last integral is convergent. The integral 

/i(dr) 
[o,i[ 1 + ^ 

is convergent, hence the inequality ([2]) is satisfied. 
Assume now that ([2]) is satisfied. Note that 

,(dr)<2/ ^<oo 

[0,1] ^[0,1] 

It follows that the integral 

1 - e"^' 



u(dr) 

[0,1] ^ 

is convergent because the integrand is bounded. It remains to consider the 
integral on [1, oo[ : 



l,oo 



[i,oo[ ''^ ^[i,oo[ ^ + r 



hence / is integrable over [0,1]. □ 
E. Proof of the CM property of the relaxation moduli in Sec. I12L 



We prove here that is CIvI for < a, 7 < 1. According to Theo- 



0,7 



rem 2.6 in Chapter 16 of j42|] it suffices to prove that F{p) := ^ pG™{p) > 

for p G C"*" and G^^{p) > for p G IR+. The second inequality follows 
from the assumption that 6 < 1. The first one follows from the fact that 
Q!,7 < 1, P £ C"*" imply that (1 + (rp)")'^ G C~ and therefore F{p) = 
—b'^ [(1 + (rp)"')^] / |1 + (tp)"!^'^ > 0. The remaining hypotheses of Theo- 
rem 2.6 are easy to verify. 

The CIvI property holds for a < 1 and 7 < 1, hence it also applies to G^^ 
and GS^. 
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